LPSC10105 
UH-511-1152-10 



O 
(N 



O 



Effective Supersymmetry at the LHC 



Howard Baertj, Sabine KramliU, Andre Lessa 1 ^, 
Sezen Sekmeni^l and Xerxes Tata 4 'u 



CN . l Dept. of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA 

2 Laboratoire de Physique Subatomique et de Cosmologie, UJF Grenoble 1, CNRS/IN2P3, 
^ ; INPG, 53 Avenue des Martyrs, F-38026 Grenoble, France 

3 Dept. of Physics, Florida State University, Tallahassee, FL 32306, USA 
4 Dept. of Physics and Astronomy, University of Hawaii, Honolulu, HI 96822, USA 
^3. 5 Dept. of Physics, University of Wisconsin, Madison, WI 53706, USA 

> 

t^- | Abstract 

We investigate the phenomenology of Effective Supersymmetry (ESUSY) models wherein 
C*~) ' electroweak gauginos and third generation scalars have masses up to about 1 TeV while 

j — . first and second generation scalars lie in the multi-TeV range. Such models ameliorate the 

SUSY flavor and CP problems via a decoupling solution, while at the same time maintain- 
ing naturalness. In our analysis, we assume independent GUT scale mass parameters for 
third and first/second generation scalars and for the Higgs scalars, in addition to m±/2, 
tan/3 and Aq, and require radiative electroweak symmetry breaking as usual. We analyse 
the parameter space which is consistent with current constraints, by means of a Markov 
Chain Monte Carlo scan. The lightest MSSM particle (LMP) is mostly, but not always the 
lightest neutralino, and moreover, the thermal relic density of the neutralino LMP is fre- 
quently very large. These models may phenomenologically be perfectly viable if the LMP 
before nucleosynthesis decays into the axino plus SM particles. Dark matter is then an 
axion/axino mixture. At the LHC, the most important production mechanisms are gluino 
production (for nil / 2 ~ 700 GeV) and third generation squark production, while SUSY 
events rich in 6-jets are the hallmark of the ESUSY scenario. We present a set of ES- 
USY benchmark points with characteristic features and discuss their LHC phenomenology. 

PACS numbers: 14.80.Ly, 12.60.Jv, 11.30.Pb 
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1 Introduction 

Particle physics models that include weak scale softly broken supersymmetry (SUSY) are es- 
pecially compelling in that they stabilize the weak scale, even in the presence of new physics 
at much higher energy scales, such as Mpi anc k or Mgut- The simplest of such models, the 
Minimal Supersymmetric Standard Model (MSSM) [TJ [2] enjoys some compelling indirect 
experimental support in that the measured values of the gauge couplings at energy scale 
Q = Mz, when run to high energies via renormalization group evolution (RGE), nearly unify at 
Q = Mqut — 2 x 10 16 GeV, as is expected in the simplest grand unified theories (GUTs). On 
the astrophysical side, SUSY theories with a conserved -R-parity offer several candidates for the 
dark matter particle, among these: the lightest neutralino, the gravitino and the axion/axino 
sup er mult iplet. 

Along with these successes, generic SUSY models also lead to new puzzles not present in 
the Standard Model (SM). Additional (usually discrete) symmetries are necessary to prevent 
a too rapid rate for proton decay. Moreover, unless weak scale soft SUSY breaking (SSB) 
parameters of the MSSM are flavor-blind pQ or "aligned" [3J, and nearly real, the model leads 
to unacceptably large rates for flavor-changing transitions and CP violating effects. We stress 
that various (sometimes, admittedly ad hoc) mechanisms have been suggested to ameliorate 
these undesired effects which, we speculate, arise because of our lack of understanding of how 
SM superpartners feel the effects of SUSY breaking. 

Another possibility to suppress unwanted flavor-changing and CP-violating effects is to 
arrange for matter sparticles to be heavy so that their effects are sufficiently suppressed [I]. 
Large matter scalar masses also serve to suppress proton decay processes in a SUSY GUT. 
This decoupling solution to the SUSY flavor problem usually requires SSB terms of order 10- 
100 TeV [5], well beyond expectations from naturalness, which favors weak scale soft terms. 
It was noted as far back as 1986 [6], and later again in Ref. [7], that the stability of the weak 
scale to radiative corrections requires only the electroweak (EW) gauginos and third generation 
sparticles - - these couple with large strength to the Higgs sector - - to have masses up to 
0(1) TeV, while gluinos and superpartners of the first two generations (whose direct couplings 
to the Higgs sector are very small, so that these enter naturalness considerations only at the 
two-loop level) could well have multi-TeV masses. Since the most stringent constraints on 
flavour- and CP-violation come from the first two generations of quarks and leptons, such a 
mass spectrum potentially alleviates the SUSY flavor and CP problems without the need for 
undue fine-tuning of parameters. This was subsequently developed into the framework referred 
to as "Effective Supersymmetry" (ESUSY), by Cohen, Kaplan and Nelson [8], who suggested 
two different realizations of the split matter generations idea by introducing a SUSY-breaking 
sector to which the first two generations couple more strongly than the third generation. As 
a result, the matter scalars of the first two generations acquire larger SUSY-breaking masses 
than third generation scalars. 

The ESUSY scenario seems in recent years to have become less favored due to two mea- 
surements. The first -- the measured (2-3)<r deviation in (g — 2)^ from its Standard Model 
(SM) expectation [S] - - seems to require smuons/muon sneutrinos in the sub- TeV range if the 
deviation is to be attributed to SUSY. The second - - the increasingly precise measurement 
of the dark matter (DM) density of the Universe - - is difficult to reconcile in ESUSY if the 



dark matter is assumed to be dominantly composed of thermal relic neutralinos left over in 
standard Big Bang cosmology. With scalars in the multi-TeV range, along with a bino-like 
neutralino, the relic density is calculated to be typically several orders of magnitude higher 
than the experimentally observed value |10j . 



fl BM h 2 = 0.1123 ±0.0035, (1) 

although there are some special parameter regions where this need not be the case. 

At the present time, it is not completely clear whether the (g — 2) M anomaly is real. The 
current discrepancy arises if one adopts the (more direct to use) hadronic vacuum polarization 
amplitude from low energy e + e~ — > hadrons data. If instead the vacuum polarization is taken 
from r lepton decay data, then the discrepancy is smaller (but recently growing!). 

In the case of neutralino DM, it is possible to have a small superpotential /i term co- 
existing with large scale masses, as occurs e.g. in the hyperbolic branch/focus point region of 
the mSUGRALJ model. Then the lightest neutralino, instead of being typically bino-like with a 
small annihilation cross section and concomitantly large relic density, becomes mixed higgsino 
dark matter with thermal relic neutralinos from the Big Bang making up the observed cold 
dark matter relic density. Other possibilities for obtaining the right relic density are resonant 
annihilation through the pseudoscalar Higgs, or co-annihilation with third-generation sfermions. 

Regions of parameter space where the neutralino relic density is too large cannot be un- 
equivocally excluded in extensions of the model. For instance, if one invokes the Peccei-Quinn- 
Weinberg-Wilczek (PQWW) solution to the strong CP problem [11], then one expects the 
presence of an axion/axino supermultiplet in SUSY theories. If the axino a is the lightest 
SUSY particle (LSP), then Z\ — > a / -f and other decay modes are allowed, which can greatly 
reduce the DM abundance far below the level expected from neutralinos. The scenario of mixed 
axion/axino DM has been examined in the SUGRA context recently in Ref. [12] in a general 
19 parameter MSSM, and appears to be at least as viable as the case with neutralino DM. 

In light of these considerations, we feel it would be fruitful to re-visit some of the phe- 
nomenological implications of ESUSY models at the beginning of the LHC era. In our analysis, 
we subsume the qualitative features of the ESUSY model: i.e. third generation, Higgs sector 
and EW gaugino masses at or below the TeV scale, with multi-TeV SSB parameters for the 
first two generations. For simplicity we also assume gaugino mass unifications While we adopt 
the qualitative picture of Effective SUSY, we do not assume the validity of the more speculative 
mechanisms that the authors of Ref. [5] suggest for the hierarchy between the SSB parameters 
of the first two generations and the corresponding parameters for the third generation and the 
gaugino sector. For our phenomenological analysis, we simply assume that such a hierarchy 
occurs for SSB parameters renormalized at a high scale that we take to be Mgut- We also 
take the less ambitious view that the SUSY flavour and CP problems are ameliorated, but 
not completely solved, by this hierarchy, and allow the first two generations of scalars to have 
masses in the 5-20 TeV range, so that some degree of universality /alignment (for the first two 



1 Minimal Supcrsymmctric model with Universal soft terms at the Gut scale and RAdiative electroweak 
symmetry breaking. 

2 While the gluino mass parameter can, in principle, be hierarchically larger, RGE effects due to a very heavy 
gluino would raise third generation squark mass parameters to high values. Gauge coupling unification also 
prefers that the gluino not be too heavy. 



generations) is still necessary to satisfy the most stringent flavour constraints. For simplicity, 
we will assume GUT scale scalar mass universality in the subspace of the first two generations, 
but allow independent SSB parameters for the third generation and the Higgs sector. This 
assumption of diagonal GUT scale scalar SSB mass squared matrices undoubtedly has an ef- 
fect on flavor physics [13] , but should have very limited impact on the implications for collider 
physics and cosmology of the ESUSY scenario — which are after all our main interest in this 
paper. 

The ESUSY model parameter space we will examine is thus given by the set of parameters 
(renormalized at the GUT scale) 

m (l,2), m (3),m Hu , m Hd , A , m 1/2 , tan/3, sign(/x) (2) 

along with the top quark mass m t = 173.1 ±1.3 GeV. Here mo(l, 2) and mo(3) are the masses of 
the first/second and of the third generation sfermions, respectively; m# u = m 2 Hu / \ \m 2 Hu \ and 

rriH d = m H d / \/\ m H d \i are the SSB mass parameters of the up- and down-type Higgs scalars; A 
is a universal trilinear coupling (relevant mostly only for the third generation); m 1 / 2 a universal 
gaugino mass parameter; and tan/3 = v u jv^. Radiative electroweak symmetry breaking can be 
used to determine /i 2 using the measured value of Mz- 

We wish to maintain the successful gauge coupling unification at Q = Mqvt, so that 
we will assume here that the MSSM is the correct effective field theory between Mqut and 
M = mo(l, 2), the scale at which first and second generation scalars decouple from the theory. 
Below M, we have ESUSY - the SM with two Higgs boson doublets together with third 
generation scalars, gauginos and higgsinos — as the effective theory. 

Models with third generation scalar and gaugino-higgsino sector at ~TeV but multi-TeV first 
and second generation mass parameters have been investigated in the past. Shortly after Cohen 
et al. [8] laid down their framework, it was noted that two loop RGE effects arising due to heavy 
first/second generation scalars act to suppress third generation scalar mass parameters (even to 
tachyonic values), and a variety of flavor and CP- violating constraints were examined IT3] . 
In Refs. [151 EH] it was shown that the inverted scalar mass hierarchy which is the hallmark 
of the ESUSYjscenario could emerge dynamically in models with Yukawa coupling unification, 
although the M : M weak ratio was found to be limited [16] • This class of models -- requiring 
t — b — T Yukawa coupling unification and SO (10) like boundary conditions with tan/3 ~ 50 - 
has been investigated in detail in Refs. [TTJ [18]. In Ref. [19], the magnitude of the M : M wea k 
hierarchy was investigated in models without Yukawa coupling unification where the matter 
scalar mass parameters are already taken to be split at the GUT scale. This study assumed 
equal GUT scale values of third generation and Higgs SSB masses, and as a result, was limited 
in scope compared to the results to be presented here. In addition, several authors have 
investigated other related aspects of supersymmetric models with heavy scalars |20j . 

The remainder of this paper is organized as follows. In Sec. [2j we explore the parameter 
space in Eq. ([2]) and map out ranges of parameters that potentially lead to ESUSY at the 
weak scale. In Sec. [3J we perform a Markov Chain Monte Carlo (MCMC) analysis to search for 
SUSY scenarios which fulfill the ESUSY conditions, subject to various experimental constraints. 
We first describe the setup of our MCMC, and then show results for posterior probability 



distributions both for input parameters that result in ESUSY, as well as for expectations for 
various sparticle masses and selected experimental observables. Moreover, we pick out some 
characteristic benchmark points for further study. In Sec. HI we first present the main particle 
production rates and decay patterns relevant for LHC phenomenology of the ESUSY scenarios, 
and then illustrate the diversity of LHC phenomena through more detailed discussion of the 
phenomenology of the benchmark points, including the possibility that the lightest MSSM 
sparticle (LMP, as distinct from LSP) may be charged or coloured. We conclude in Sec. [5] with 
a summary of our results. 

2 Viable spectra and parameter space of the ESUSY 
model 

Constrained supersymmetric models such as the mSUGRA model have, for a fixed gluino mass, 
an upper limit on how massive scalars can be. For a model with universal scalar mass soft 
terms at the GUT scale, given by m , with m.1/2 fixed, the weak scale value of n 2 diminishes 
as m increases. Ultimately, when m is large enough, |/x| becomes comparable to the weak 
scale value of the bino mass Mi, and one gains neutralinos of mixed higgsino-bino composition 
which make a good candidate for thermal WIMPs. This is the so-called hyperbolic branch/focus 
point (HB/FP) region of the SUGRA parameter space [21]. As m increases even further, |/x| 
decreases even more and the Zi becomes nearly pure higgsino, and the higgsino-like lightest 
chargino may drop below limits from LEP2 searches. For yet higher values of mo, /x 2 becomes 
negative, signaling an inappropriate electroweak symmetry breaking (EWSB) pattern. Thus, 
in the mSUGRA model, for a given value of mi/2, mo can become no larger than a few TeV. 

If we proceed to models with non-universal generations, taking mo(l,2) as independent 
of m (3), with m (3) = m# u = m# d = m# as in Ref. [19], we expect, using 1-loop RGEs, 
by the same reasoning, an upper bound on m (3). The point is that the first and second 
generations essentially decouple from the RG evolution of the Higgs SSB parameters (which 
enter the electroweak potential minimization conditions, see Ref. [22]). At two loop order, 
however, very large scalar masses do affect the running of the other scalar masses and in fact 
lead to their suppression at the weak scale [TU [19]. As m (3) decreases (or m (l,2) increases), 
so m (3) <C m (l, 2), one of the third generation scalars becomes the LSP, turning tachyonic for 
even smaller (larger) m (3) (mo(l, 2)) values. Therefore, with t«q(3) = m 2 H at the GUT scale, 
there are both an upper (from /i 2 > 0) and a lower (from m?r _ > 0) bound on m (3). In Fig. [U 
we show regions of allowed parameter space in the m# versus m (3) plane for mo(l, 2) = 5, 10 
and 20 TeV, with a) m.1/2 = 300 GeV and b) m.1/2 = 800 GeV, where the spectrum calculations 
are performed using Isajet7.80 [23]. The green dashed line is where tuh = mo (3). The region 
to the right of the blue lines is forbidden due to lack of EWSB, while the region below the red 
lines is forbidden because a third generation scalar becomes tachyonic. The region to the left 
and above the contours yields viable spectra. 

For low mi/2, we see from Fig. [T]a) that if we restrict our attention to m# = m (3) (green 
dashed line) only the interval 1 TeV < m (3) = m H < 4 TeV is allowed for m (l,2) = 5 
TeV, while there are no allowed values for mo(l,2) = 10 or 20 TeV, since for these cases the 
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Figure 1: Regions of parameter space leading to viable sparticle mass spectra in the tjih vs. m (3) 
plane for m (l,2) = 5, 10 and 20 TeV. We fix A = 0, tan/3 = 10 and \i > and choose 
m i/2 = 300 GeV in the left frame andmi/ 2 = 800 GeV in the right frame. The region below the 
red lines is excluded because of tachyonic third generation masses or, for large negative values 
of m H because m 2 A < 0, while the region to the right of the blue lines does not allow EWSB. 



The green dashed line has m (3) 
lowest allowed values. 



m H . The ESUSY spectra results from taking m (3) to its 



mo (3) = rax line lies entirely in the excluded region. However, through 1-loop RGE running 
effects, an increase in mi/2 increases rriq. Hence larger values of m (l,2) become viable if 
mi/2 is large enough. This is seen in Fig. [Tf)), where the green dashed line (again representing 
mo(3) = ran) shows that, for mi/ 2 = 800 GeV, the range of allowed m (3) = ran values is larger 
for mo(l, 2) = 5 TeV, where the tachyonic lower bound is gone, so < mo(3) = nin ^ 5 TeV. 
Besides, if a neutralino LSP is required, we have 0.25 < mo(3) = m# < 5 TeV. For the higher 
mo(l, 2) cases we see that the increase in mi/2 provides an allowed region for mo(l, 2) = 10 TeV, 
but still is insufficient for the mo(l,2) = 20 TeV case. It is for this reason that the solutions 
found in Ref. [IH] with very large m (l,2) values, as needed to suppress FCNCs, also required 
rather high m^ values, and consequently very heavy gluinos: quite beyond the reach of LHC. 
Since there is no compelling theoretical argument to link the Higgs and third generation 
mass scales, it seems reasonable to adopt independent values of mo(3) and mn u = m H d = mil 
(NUHM1) or rriH u ^ m H d (NUHM2). In this case, for a large value of mo(3), we can take 
m 2 H <C rnl(3), and thus give m H a head start on running towards the necessary negative squared 
values which are needed for successful EWSB (recall that at the weak scale, /i 2 ~ — m 2 H for even 
modest values of tan/3). This is the well-known feature of NUHM1 models wherein increasing 
rriH results in a decrease of the weak scale value of |//| [21]. Note that large negative values of 
ran are excluded because there m 2 A < 0, signaling that EWSB is not correctly obtained. 



m (l,2) = 10 TeV, m = 0.5 TeV, n^ = 0, A Q = 1 TeV, tar(3 = 10, |l > 0, rn = 173.1 GeV 



1500 







1000 



mJ3)(GeV) 



Figure 2: Contours of fi, gluino, stop and Higgs (h) boson masses in the m (3) vs. mi/2 plane 
for m (l,2) = 10 TeV, m Hu = 0.5 TeV, m Hd = 0, A = 1 TeV, tan/3 = 10 and fi > 0. 
Also shown are regions where various scalars are the lightest MSSM particles (LMPs). In the 
almost impossible to see narrow green band whose thickness is ~ 10 GeV that hugs the upper 
boundary of the stop LMP region, the thermal relic abundance Q^ h 2 < 0.11. The gray shaded 
region is excluded because either we have tachyonic sparticles or the chargino mass is below the 
corresponding LEP2 limit. 



If we adopt a low value of |m#| ~ 1 TeV, and then move from high to low m (3) values, 
the third generation scalars decrease in mass. The region for ESUSY with third generation 
scalars at ~TeV is at mo (3) values just above the tachyonic third generation region. Thus, the 
strategy for gaining viable ESUSY spectra is to 1.) adopt a large value of m (l,2) ~ 5 — 20 
TeV, then 2.) adopt a low value of mg ~ — 2 TeV, and finally, 3.) starting at a several TeV 
value of mo (3), decrease its value until third generation scalars dip below the TeV region. If 
we want low |/x| values as well, then increase m# as close to the EWSB boundary (on the right 
of the figure) as desired. 

It is instructive to portray ESUSY in various parameter planes. Figure [2] shows the allowed 
region in the m (3) vs. m 1 / 2 plane with m (l, 2) = 10 TeV, m Hu = 0.5 TeV, m Hd = TeV (along 
with A = 1 TeV, tan/3 = 10 and /i > 0). The left-side gray shaded region gives tachyonic stop 
or sbottom masses, while in the lower-right gray region the chargino mass is below its LEP2 
limit. The white region gives viable supersymmetric mass spectra. We plot contours of m~ g 
(gray dashed), m^ (blue solid), m^ (dashed cyan) and \i (magenta dot-dashed). The region to 
the lower left of the cyan dashed contour has m^ < 114 GeV. The slim region adjacent to but 
right of the tachyonic stop region gives viable ESUSY spectra with top squark masses mi < 1 

TeV. There is an almost-impossible-to-see green region of thickness ~ 6 GeV in mi/2 hugging 
the boundary of the tachyonic stop region where the thermal neutralino abundance is in accord 



with measurement: Q?? h 2 < 0.11. This is, in fact, the top squark co-annihilation region [23] for 
ESUSY. Along the tachyonic boundary is another narrow region where the t\ can become the 
lightest MSSM particle (LMP), though not necessarily the LSP if a non-MSSM i?-odd sparticle 
(such as the axino) is lighter. We also show the regions where 61, f\ or v t l is the LMPj; 



In Fig. |3]a), we fix instead the values of m (l, 2), m (3) and m 1( / 2 and plot contours of 
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/j, and VfiA in the rriH u vs. rn,H d plane. As discussed above, for intermediate positive ranges of 
m,H d in the plot, where mn u ~ ™>H d , as m H u increases, the value of \x decreases until /i 2 < and 
EWSB is no longer realized. This can be seen by the /1 contours in this part of the figure where 
we also see a significant green region where the thermal neutralino relic density is consistent 
with its observed value. In the lower portion (low m# d ) of the grey area on the right and at 
the bottom in Fig. [SJa), we do not obtain radiative EWSB because rn? h or m 2 A become negative. 
The green region in this vicinity is where the neutralino relic density is compatible because of 
higgs resonance annihilation [SB]- Finally, for large enough positive values of m# u and mH d , t\ 
becomes the LMP (shaded blue region), with m~ < above that (upper portion of the grey 
area on the right) and co-annihilation with stops leads to a compatible neutralino relic density. 
Finally, we present in Fig. [3]&) the allowed region in the m# u vs. m (3) plane with mo(l, 2) = 10 
TeV, mi/2 = 0.5 TeV and m Hd = 1 TeV. We see that the ESUSY scenario is realized at low 
m (3) and m Hu < m (3). 

3 Markov Chain Monte Carlo analysis of ESUSY 

In this section, we perform an MCMC analysis to map out the regions of ESUSY parameter 
space that are consistent with current experimental constraints from collider and flavour physics. 
The characteristics of ESUSY are enforced by means of theoretical priors. We first briefly 
describe our set up and list the likelihood functions and model priors that we use, and then 
show our results for the posterior probability distributions both for input parameters as well 
as various sparticle masses and selected experimental observables. In Sec. 13.31 we present some 
benchmark points for more detailed study. 

3.1 Setup of the MCMC scans 

The setup and procedure of our MCMC closely follow [22], so we do not repeat these here 
but instead just give the key data of the analysis. For details on the MCMC method and 
Bayesian analysis, see e.g. [281 EHl EQ]. In addition to Isajet7.80 [23], which we use for 
the spectrum calculation, we use SuperIso2.7 [ST] for computing B-physics observables and 
micrOMEGAs2.4 1 32 1 for the calculation of the neutralino LMP abundance. 



3 With a common value of too (3), the v T LMP may be surprising since at the one- loop level Yukawa inter- 
actions drive m- to smaller values than m% . In this case, however, two-loop effects due to heavy scalars are 
large, and for smaller values of mii 2 can lead to rrif L < Wf B and a sneutrino LMP because of the D-term. Since 
one-loop effects increase with mw2, ultimately the usual situation with stau LMP is obtained. Incidently, the 
compensation between the one and two-loop contributions implies that for certain parameters we can have a 
near equality of the left- and right-stau mass parameters, and so large stau mixing even though the tau Yukawa 
coupling is not particularly large. 
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Figure 3: The uppper frame a) shows contours of \x, m^ and m^ in the vrtn u vs. mH d plane 
for mo(l,2) = 10 TeV, m (3) = 3 TeV, m 1/2 = 0.5 TeV, A = 0, tan/3 = 10 and fi > 0. In 
the green band, the thermal relic abundance Q^ h 2 < 0.11 and in the shaded blue region the t\ 
is the lightest MSSM particle. The gray shaded region is excluded as explained in the text. The 
lower frame b), shows contours of fj, and m^ in the m Hu vs. m (3) plane for m (l,2) = 10 
TeV, mi/2 = 0.5 TeV, m Hd = 1 TeV, A = 0, tan/3 = 10 and // > 0. 



For each point P in the scan of ESUSY parameter space, we first compute individual 
likelihoods L(Oi) for each experimental observable 0{. Then the overall likelihood is given as 
the product of these individual likelihoods as Lp = rjL(Oj). The observables we use are listed 
in Table [1] along with the parameters of the corresponding likelihood functions used. The forms 
of these likelihood functions are given by, 



Lx(x,x ,a a 



exp[ 



[x 



x )/cr x 



L 2 (x,x ,a x ) =exp 






(3) 



Observable 


Limit 


Likelihood function 


Ref. 


BF(b -> S7) 


(3.52 ±0.34) x 1(T 4 


L 2 (x,3.52x 10~ 4 ,0.34x 10" 4 ) 


[Mi Eg 


BF(£? S -> //+//-) 


< 5.8 x 10" 8 


Li(x,5.8 x 10" 8 ,-5.8 x 10" 10 ) 


[36] 


#(£„ -> rz/ T ) 


1.28 ±0.38 


L 2 (x, 1.28, 0.38) 


[34] 


m t 


173.1 ±1.3 


L 2 (x, 173.1, 1.3) 


[37] 


m h 


> 114.5 


1 or exp(-x 2 (m h ,)/2) 


[33] 


SUSY mass limits 


LEP limits 


1 or 10" 9 


[38] 



Table 1: Observables used in the likelihood calculation. 



for observables for which there is only an upper or lower bound, and for observables for which 
a measurement is available, respectively. For the LEP limit on the Higgs mass, we use [33] 



X 2 (m h ) 



(m h 



m 



limit \ 2 



1.1 GeV) 2 ± (1.5 GeV)< 



(4) 



with m 



limit 
h 



115 GeV. The likelihood L{rrih) is then given by Lim^) = e x t" 1 '*)/ 2 for m^ < 
115 GeV and L{m h ) = 1 for m h > 115 GeV. 

In order to favour (sub)TeV-scale electroweak gauginos, higgsinos, and third generation 
sfermions, we multiply the likelihood obtained from the experimental constraints with the 
following model prior: 



L MsS ( m x) 



1 



exp((m x - M eff )/170 GeV) 



for each X = xt, X2 > ^i> ^2, b±. Choosing, for instance, M e s = 1.5 TeV, we get L 



1\L. 



(5) 



0.95, 
0.95, 



0.5, 0.05 for mx = 1, 1.5 and 2 TeV respectively; for M e g = 1 TeV, we get Lm cS 
0.5, 0.05 for mx = 0.5, 1, 1.5 TeV respectively. We have run chains for different M e g and 
checked that the effect is indeed only to vary the upper bound on the effective SUSY masses; 
the qualitative features of the parameter space remain unchanged. 

Finally, regarding the model parameters, we allow myz — [0, 2] TeV, m (3) = [0, 10] TeV, 
m Hud = [-10, 10] TeV (with m Hud unrelated), A = [-40, 40] TeV, and tan/3 = [2, 60], taking 
flat prior probability distributions for all these parameters. In addition we let the top mass 
vary within m t = 173.1 ± 1.3 GeV [39] with a Gaussian distribution. For m (l,2), we show 
results for two different approaches: 

1. we let mo(l, 2) vary from 5 to 20 TeV with a flat prior probability distribution (L^j = 1), 
or 



2. we let m (l, 2) vary without limits but apply a model prior of 

1 



L^(m (l,2)) 



1 ± exp((10 - m (l, 2))/1.7 TeV) 
in order to favour m (l, 2) values beyond 10 TeV. 



(6) 
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The total likelihood of a point is then taken to be L tot = Lp x Lm cS x Lj^. 

Before turning to the results, a comment is in order regarding the nature of the LMP. If 
we do not impose any dark matter requirement on the LMP, the large majority of points have 
a neutralino LMP with a far too high relic density of up to Q x h 2 ~ 10 3 . As noted in the 
Introduction, this can be in agreement with cosmological observations if the "true" LSP is 
actually an axino, i.e. in the case of mixed axion/axino dark matter [121 ED]- Moreover, about 
12% (4%) of the points accepted by the chains with M eS = 1 (1.5) TeV have a stop, sbottom 
or stau LMP. Again these points may be viable if the true LSP is the axino. However, as we 
will discuss below, the phenomenology in these charged LMP cases is quite different from the 
neutralino LMP case: in particular, if the stop, sbottom or stau LMP does not decay promptly 
but gives a charged track. In what follows, we will confine our MCMC analysis to the case of 
a neutralino LMP, which indeed occurs most frequently. We will, however, comment on the 
phenomenology of the colored or charged LMP case at the end of the next section. 

3.2 MCMC results with neutralino LMP 

Figure H] shows the posterior probability distributions of the input parameters of the ESUSY 
model from MCMC scans requiring a neutralino LMP. Here we have used M c g = 1 TeV in 
Eq. (JSJ). Unseen dimensions are marginalized over. The figure compares the two priors for 
mo(l, 2): The thin black lines are for case 1. where m (l, 2) = 5—20 TeV with a flat distribution, 
while the thick red lines are for case 2. which has no limits on m (l,2) but the choice of the 
prior in Eq. §6§ favours higher values of m (l, 2). We see that, as anticipated, we find solutions 
where third generation sfermion and gaugino SSB parameters are typically 1-2 TeV, while the 
first/second generation SSB parameters are favoured to lie beyond 10 TeV. Moderate values of 
tan/3 are favoured. 

As can be seen, the precise condition on mo(l, 2) hardly influences the probability distri- 
butions of the other parameters. A possible exception is mi/2, which features a double-peak 
distribution. This comes from the fact that at small mi/2, the parameter space volume is con- 
strained in the directions of the scalar mass parameters, while there is more space in the tan /3 
direction. For large mi/2, on the other hand, tan/3 is very much constrained by BR(6 — > 57), 
while there is more volume in the large mo (1,2) directions. The effect is more pronounced 
for case 1. which does not disfavour smaller values for m (l,2). Moreover, we see that larger 
mo(l,2) prefers somewhat larger m^Jj It is also interesting to note that very high mo(l,2) 
around 15-20 TeV suffers from a low probability. This is because as m (l, 2) increases, we are 
forced to increasing larger values of mo(3) (see Fig. [1]). For values of mo(3) in the multi-TeV 
range, third generation SSB parameters at ~TeV can occur only when the large two loop RGE 
contributions rather precisely cancel the naturally multi-TeV contribution that we start with - 
too little cancellations leave large SSB parameters, while too much cancellation leads to tachy- 
onic masses. As a result, the region of parameter space with a light third generation rapidly 
shrinks when mo(l,2) exceeds ~ 15-20 TeV. 



4 The shape of the posterior probability distribution of m.1/2, and tan/3 are likely the most sensitive to the 
b — > s"/ constraint, and hence to our assumption that there are no off-diagonal pieces in the GUT scale squark 
mass matrices. 
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Figure 4: Posterior probability distributions of the input parameters of the ESUSY model from 
MCMC scans with M e s = 1 TeV, comparing the two priors for mo(l, 2). We take fj, > 0, and 
assume that the neutralino is the LMP. The thin black lines are for case 1. where mo(l,2) = 
5 — 20 TeV with a flat distribution, while the thick red lines are for case 2. which has no limit 
on m (l, 2) but the prior Eq. |6j) disfavours low values for this parameter. The last panel shows 



the combination m 2 



,(1)2 



,(1)2 



The posterior probability distributions of several relevant masses are shown in Fig. [5j We 
infer from the EW gaugino spectrum (and also the //) distribution that most of the time the 
lightest neutralino is bino-like. 

Figure [6] displays the probability distributions of £>-decay observables, Aa M — (g — 2)^ USY , 

the higgsino fraction f H 



+ t>2 (in the notation of Ref. [22]) of the Z 1; and the relic 
abundance Qh 2 of the Z\. It is clear that if the muon anomalous moment is confirmed, it 
cannot arise within the ESUSY framework. We note that there is non-negligible fraction of 
the parameter space with a higgsino-like LMP. We also see that the neutralino relic density 
peaks around Qh 2 ~ 1 — 10 and goes up to Vth 2 ~ 10 4 . Nevertheless, solutions with low values 
of Vth 2 ~ 0.1, where the relic density of thermal neutralinos does not yield a universe that 
is too short-lived, also have non-negligible probability. They occur because of Higgs funnel 
annihilation, coannihilation with light stops and/or sbottoms (or, less likely, staus), or because 
of a large LMP higgsino component. 



3.3 Effective SUSY benchmark points 

To illustrate the types of SUSY spectra and the diverse phenomenology that can result in 
ESUSY models, we list five benchmark points in Table [2J The first three points we list have 
a neutralino LMP, while points ES4 and ES5 have a sbottom and a stau LMP, respectively. 
Point ESI adopts an mo(l,2) = 10 TeV, which is the most favored value in the probability 
distributions resulting from the MCMC. On the other hand, ESI has a value of m~ g = 524 GeV, 
not far from the ESUSY minimum of nig ~ 460 Gev. The chargino W\ and the neutralino Z 2 
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Figure 5: Same as Fig. [£] but for various sparticle and Higgs masses, and [i. 



have masses of just 139 GeV, and the \i parameter is 857 GeV. All third generation squarks 
except b 2 , as well as all EW gauginos are ~1 TeV scale, while the two staus are over 2 TeV. 

In point ES2, mo(l, 2) is increased to 12 TeV. ES2 has much higher value of m§ ~ 2.4 TeV 
with concomitantly heavier W\ and Zi j2 , although all EW gauginos along with t\, i 2 and b\ are 
below 1 TeV. The b 2 and staus all have masses ~ 1.3 — 1.4 TeV range. For this point, t\ decays 
exclusively via the three-body mode t\ — > bWZi, while other third generation quarks decay via 
two-body decays. 

Point ES3, with a value of mo(l,2) ~ 10 TeV, again has all third generation sfermions, 
as well as all charginos, neutralinos and Higgs bosons essentially at or well below 1 TeV. The 
gluino mass is rather high: ~ 2.1 TeV. One characteristic feature of this point is that, because of 
the near degeneracy of £1 and Z 1; not only the two-body decays but even the three body decay 
£1 — > bWZi is kinematically forbidden. In this case, £1 will decay via ti — > cZ\ or t\ — \ bZiff, 
though of course the small available phase space will favour the former decay [H] . 

Point ES4 has a high mo(l,2) = 20 TeV, a heavy gluino with a mass of 2.5 TeV, and 
correspondingly heavy EW gauginos. However, t\ and b\ are very light with masses of 327 GeV 
and 291 GeV, while m(i 2 ) = 793 GeV. In this case, b x is the LMP. 

Finally, point ES5 has been chosen to illustrate that staus can also be very light. Here, f\ is 
the LMP, with a mass of just 289 GeV, while m(f 2 ) = 380 GeV. The lighter stop and the EW 
gauginos are in the sub- TeV range, while other third generation sfermions are 1.1-1.2 TeV. 

Two of the ESUSY points with a neutralino as LMP (ESI and ES2) yield a thermal abun- 
dance of neutralinos far in excess of WMAP measured value for the cold DM relic density, while 
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1.1 
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1 
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-0.6 


1 





m H u [TeV] 


2 


1.2 


0.5 


1 





mi/a [TeV] 


0.16 


1 


0.85 


1 


0.8 


^o [TeV] 





0.9 


1 








P 


857 


902 


925 


2329 


855 


m- g 


524 


2446 


2103 


2526 


1941 


rriu L [TeV] 


10.0 


12.1 


10.1 


20.1 


6.2 


rriu R [TeV] 


10.0 


12.2 


10.1 


20.1 


6.2 


mg L [TeV] 


10.0 


12.0 


10.0 


20.0 


6.0 


m e - s [TeV] 


10.0 


12.0 


10.0 


20.0 


6.0 


"i*i 


646 


608 


398 


327 


867 


m t 2 


1049 


948 


770 


793 


1147 


m ' bl 


1039 


830 


586 


292 


1083 


m h 


1711 


1313 


958 


1885 


1176.3 


m fl 


2269 


1341 


944 


2956 


289 


m f2 


2299 


1388 


1008 


3152 


381 


m Wl 


139 


815 


708 


881 


658 


m z 2 


139 


815 


708 


878 


657 


m Zl 


69 


441 


372 


452 


347 


m A 


1022 


450 


398 


2042 


875 


m h 


110.7 


118.3 


117.3 


119.9 


117.8 


tt^h 2 


320 


0.789 


0.036 






cr[fb] 


23.2 x 10 3 


157.8 


1618.0 


12.3 x 10 3 


51.4 


99 


62.8% 


0.02% 


0.01% 


- 


1.5% 


EW — ino pairs 


36.6% 


5.1% 


0.8% 


0.06% 


35.3% 


slep. pairs 


- 


0.02% 


0.01% 


- 


24.4% 


% 


0.4% 


78.3% 


87.6% 


28.2% 


26.5% 


hh 


0.03% 


4.5% 


1.8% 


0.3% 


3.8% 


hh 


0.02% 


11.5% 


9.1% 


71.4% 


4.9% 


b 2 b 2 


- 


0.4% 


0.4% 


- 


2.9% 



Table 2: Masses and parameters in GeV units (unless otherwise noted) for five case study points 
in ESUSY using Isajet 7.80 with m t = 173.1 GeV, tan/3 = 10 and /i > 0. The LMP mass is 
denoted with bold numbers. We also list the total tree level sparticle production cross section 
in fb at the LHC with a center of mass energy of 14 TeV. 
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Figure 6: Same as Fig. [7] but for a variety of low energy observables, the higgsino content of 
the lightest neutralino, and the thermal neutralino relic density. 



point ES3 is in the top-squark co-annihilation region and gives a thermal neutralino abundance 
of Qz h 2 = 0.036. As mentioned previously, the points ESI and ES2 can still be cosmologically 
viable if we invoke the PQWW solution to the strong CP problem, which necessitates the in- 
troduction of an axion/axino supermultiplet into the theory. In this case, the axino a, and not 
the neutralino, can be the LSP. If the a is the LSP, then the neutralinos will decay via Z\ — \ d'-f 
with lifetime smaller than ~ 0.1-1 s, i.e. before the onset of Big Bang nucleosynthesis (BBN), 
for a PQWW scale /„ ~ 10 12 GeV [12]. Since every neutralino decays to an axino, the would 
be thermal relic density of neutralinos is reduced by a factor rria/m^ , and for small enough 
rria would be compatible with the relic density measurement. Thermal production of axinos 
in the early universe can still proceed, but the relic abundance of axinos then depends on the 
re-heat temperature Tr after inflation. The dark matter will then consist of an axino/axion 
admixture [12], with relic axions being produced via the vacuum mis-alignment mechanism. 

In the cases where t\, b\ or f\ are LMPs, such scenarios are again allowed if, as before, we 
assume that the axino is the LSP. In these cases the sfermions dominant decay is f\ — > fa, 
via loop-mediated processes; the competing three-body decays, f\ — > raj [13] or q — > qag [H] 
are argued to have a small branching fraction. The loop calculation is complicated by the fact 
that the non-renormalizable axino-bino-photon (and the axino-bino-Z and axino-gluino-gluon) 
coupling enters the calculation. Nevertheless, these authors find that for m f = 100 GeV, the 
stau LMP lifetime ranges from about 0.01s to 10 h, for the PQWW scale f a in the range 
5 x 10 9 -5 x 10 12 GeV, and scales inversely as the stau mass. The corresponding squark lifetime 
is shorter: about 2 x 10~ 6 s for 500 GeV squarks and 1 TeV gluinos with f a = 10 11 GeV. If the 
LMP lifetime indeed exceeds a few seconds, this could disrupt the successful predictions of Big 
Bang nucleosynthesis [IS], though it appears that if f a < 10 12 GeV, a stau LMP is relatively 
safe [16] . For a general discussion of charged LMPs, see Ref . [17] . 



14 



pp@14TeV ° L V) 



g(tt + bb) | 




pp@14TeV - L0 (fb) 




1500 2000 2500 3000 3500 4000 4500 5000 
mJ3) (GeV) 



1500 2000 2500 30O0 3500 4000 4500 5000 
ml 3) (GeV) 



Figure 7: a) Gluino pair production and b) third generation squark pair production cross sections 
at LHC with y/s = 14 TeV in the mo(3) vs. min plane for m (l, 2) = 20 TeV, m^ j - TeV 
and Aq = 0, tan/3 = 10 and fi > 0. 
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4 Phenomenology of Effective SUSY at the LHC 

How will effective SUSY manifest itself at the LHC? To answer this, we first show in Fig. [7] 
total gluino and third generation squark pair production cross sections at a 14 TeV pp collider 
in the mo(3) vs. mi/2 plane displayed in Fig. [2j but now with mo(l,2) = 20 TeV. In frame 
a), the pp — > ggX reaction (where X denotes assorted hadronic debris) reaches to over 10 4 fb 
for gluino masses as low as m g ~ 500 GeV. The gluino pair cross section drops continuously as 
mi/2, or alternatively m§, increases. 

The ESUSY region with third generation squarks having masses ~TeV lies adjacent to 
the boundary of the gray tachyonic region. In frame b), we see the pp — > t{ti plus b^bi (for 
i — 1 — 2) summed cross sections. These are typically in the 1-100 fb range all along the 
tachyonic boundary, with portions reaching cross sections as high as 10 3 fb. Thus, the low mi/2 
region will consist of mainly gluino pair production, where the main influence of light third 
generation squarks will be upon the gluino branching fractions. As we move to higher m\ii 
values, the gluino pair cross section will diminish, and the total SUSY production cross section 
will be dominated by pair production of third generation squarks. 

Unless the gluino is very heavy we would expect that, especially after selection cuts, LHC 
phenomenology would largely be determined by gluino pair production. SUSY event topologies 
would thus be sensitive to gluino decay patterns. In Fig. [HI we show the gluino branching 
fractions for the MCMC scan points with m- g < 1 TeV. Of course, when gluino two-body 
decays are kinematically accessible, these would dominate. We see from the figure that these 
are accessible for relatively few points when m g < 500 GeV. We will return to this later when 
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we discuss how one might distinguish the ESUSY model with light gluinos from the model with 
t-b-r Yukawa unification where m g is bounded from above [17]. For the bulk of the scanned 
points, the gluino decays via three body modes into top and bottom quarks plus EW gauginos; 
corresponding decays to the first two generations are suppressed because these squarks are very 
heavy. Bottom-jet tagging will clearly provide an effective way for enriching the ESUSY event 
sample at the LHC [4*8] 09]. As in many SUSY models, direct decays to the neutralino LMP 
never have a very large branching fraction. The fact that gluino decays have large branching 
fractions to W±, Z 2 and top quarks implies that the ESUSY gluino event sample will include 
multi-lepton events from gluino decay cascades. Finally, we see from the last frame that while 
the branching fraction for the radiative gluino decays [50] g — > gZi is usually small, it can reach 
40-50%. We have checked that the neutralino in question is mostly Z 3j4: . This occurs when we 
have large splittings among third generation squarks with stops being lighter than sbottoms, 
when a light neutralino has a significant up-higgsino component and the decay g —¥ tiZi is 
kinematically forbidden. The branching fraction for decays to light quarks is very small for 
mo(l,2) > 5 TeV, although larger values are possible with the second L^ prior (case 2. in 
Sec. 13. ip . which disfavours lower values of mo(l, 2) but does not exclude them. 

4.1 Light gluinos with Z x as LMP (ESI) 

The benchmark point ESI features a rather light gluino with mass m g = 524 GeV, corre- 
spondingly light charginos and neutralinos, along with three sub- TeV third generation squarks. 
From Table [2J we see that total sparticle production cross section, summed over all reactions, 
is 23.2 pb. Of this, 62.8% comes from pp — > ggX production, while 36.6% comes from EW 
gaugino production (sum of all chargino and neutralino production cross sections). Only a tiny 
fraction of the production cross section comes from third generation squarks, so the rate for 
iiii and bibi production (for i — 1, 2) is very small. 

The gluino decays via three body modes as g — > bbZ 2 (37%), g — > ttZ\ (19%) and g — > 
tbWi + c.c (41%). Thus, each gluino pair producton event will contain at least four 6-jets, plus 
other jets, isolated leptons (from top quarks and EW gauginos) and E^ ms . The W\ decays via 
three-body modes W\ —¥ Ziff in accord with W* propagator dominance (e.g. W\ — > Z\ev e 
at 11%). The Z 2 decays to bbZ\ 20% of the time, while leptonic decays such as Z 2 — >■ Z\e + e~ 
occur at the 3% level. In this case, the b\ is quite light and dominantly bi,, which enhances the 
Z 2 decay to b quarks at the expense of first and second generation fermions. Clean trilepton 
signals from W\Z 2 production may also be observable |51| . 

The ESI scenario will have much the same LHC phenomenology as the Yukawa-unified 
SUSY with a 500 GeV gluino [49], which is also dominated by gluino pair production followed 
by three-body gluino decays into mainly 6-quarks. A natural question to ask is whether it is 
possible to distinguish these scenarios at the LHC. One difference is that Yukawa-unified SUSY 
requires a very large A Q parameter with A ~ — 2m 16 , where miQ is the matter GUT scale 
scalar mass parameter, which lies typically in the range 8-20 TeV. This large value of A feeds 
into gaugino mass evolution at two-loops in the RGEs, and suppresses m g with respect to m^ 

and m z 2 ,m- 

To illustrate this, we show in Fig. [9] a scatter plot of m g versus m^ — m^ obtained by 
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Figure 8: Results for the branching fractions for gluino decays from MCMC scan points with 
nig < 1 TeV versus the gluino mass. The upper two frames show the branching fractions 
for three-body decays of gluinos to third generation quarks plus EW gauginos. The bottom-left 
frame shows the branching fraction for two-body decays to third generation quarks and squarks, 
while the bottom-right frame illustrates that the radiative decays g — > gZi can sometimes be 
substantial. All results are for the L^ = 1 case with 5 TeV < mo(l) 2) < 20 TeV. 

scanning the parameter space of (a) the ESUSY model (red pluses) and (b) the Yukawa-unified 
(YU) model with "just-so" GUT scale Higgs soft mass splitting (HS, dark blue stars) and the 
YU model with full D-term splitting, right-hand neutrinos and 3rd generation scalar splitting 
(the DR3 model [52], light blue squares). Here, we require in YU models that R, the ratio of 
the largest to the smallest of the GUT scale Yukawa couplings, is smaller than 1.05. While 
there is significant overlap of the two classes of Yukawa-unified models, the ESUSY model lies 
on a distinct band. We expect that the neutralino mass difference edge will be rather well 
determined by the m(£ + £~) distribution. Then, even a crude measure of the gluino mass at the 
10-20% level via M e ff [53] or total cross section [51] would suffice for this distinction. Also, in 
Ref . [19] , it is shown that the mass difference m g — m^ might be extracted via the kinematic 
edge in the m{bb) distribution^ If m g can be determined even more precisely using mT2 [55] , 



We have also checked that a scatter plot of ro<j 



also serves to separate the ESUSY 



model from the Yukawa-unified SUSY scenarios, but do not show it here for brevity. 
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Figure 9: Scatter plot of the gluino mass versus the neutralino mass difference m^ — m^ that 
can be obtained from the determination of the dilepton mass edge at the LHC. We show results 
for the ESUSY model (red), and two different models with t-b-r Yukawa unification (medium 
and light blue) discussed in the text. 

its generalizations (56], or other methods that have recently been suggested, the distinction 
between the models would be even sharper. 

4.2 Light top and bottom squarks, heavy gluinos and Z\ as LMP 

(ES2 & ES3) 

For point ES2, since myz is large, the value of m g is also large, and so gg will not be produced 
at an appreciable rate at LHC. Instead, the total SUSY cross section of o~susy = 157.8 fb 
is dominated by 78.3% tit\ production, and 11.5% b\b\ production. A few percent of t 2 ^2 
also contribute. The t\ — > bWZ\ decay occurs at nearly 100% branching fraction, so the 
t\t\ production results in a bbW + W~ + E^- 1SS final state. This should be separable from ti 
background if large E^ ms or M e ff is required. Another possible background is Zti production. 
Because gluinos are heavy, SUSY contamination to the signal [57] is not an issue. The b\ — > Wt\ 
decays at nearly 100% branching fraction, so this component of the production cross section 
will result in a bbW+W'W+W- + E™ iss final state. 

The case of ES3 has four sub-TeV third generation squarks plus tau-sleptons near a TeV 
or below. The total production cross section of asusY = 1618 fb comes from 87.6% i\ti 
production, 9.1% 61&1 production plus a few percent of heavier top and bottom squark pairs. 
The visible decay products from t\ decay will be soft since there is only a 26 GeV mass gap 
between the ti and the Z\\ the visible decay products of t\ may not reliably detectable unless 
they are boosted. The b\ —¥ W~t\ at nearly 100% branching fraction, so from b\b\ production 
we expect a W + W~ + E™ ss final state, accompanied by soft charm jets or other soft hadronic 



18 



debris. The possible background comes from W + W~Z production. 

Before turning to the next benchmark case we remark that, although we have not examined 
an example of such a benchmark point, it has been shown [58J that if the gluinos and stops 
both have sub-TeV masses and that the decays g — > tti — > tbWi and/or g — > bb\ — > btW\ are 
kinematically accessible and occur with large branching fractions, a partial reconstruction of 
the event using techniques to tag the secondary top quark may be possible. Moreover, in [59] 
it has been shown that if g — > tt\ — > tcZ\ dominates, the signature of 2 b-jets plus 2 same-sign 
leptons plus additional jets and missing energy is an excellent discovery channel for gluino 
masses up to about 900 GeV. 

4.3 Bottom or top squarks as LMP: the case of heavy quasi-stable 
colored particles (ES4) 

Benchmark point ES4 illustrates an ESUSY model with t\ and b\ lighter than Z\, and in fact 
the bottom squark bi is the LMP. This is viable if the axino is the true LSP so that b\ — > ba, 
and the dark matter consists of an axion/axino admixture. For case ES4, with a PQ breaking 
scale of order f a ~ 10 11 GeV, we find a bottom squark lifetime of ~ 10~ 6 sec, i.e. the &i is stable 
as far as collider searches go, but it decays well before the start of Big Bang nucleosynthesis. 
Previous searches for quasi-stable squark hadrons Q at ALEPH find that raq > 95 (92) GeV 
for up- type (down- type) squark hadrons [60J. Searches by CDF evidently require nig ~ 241.5 
GeV [Ml E2]. At the LHC, for case ES4, sparticles will be produced with total cross section 
(Jsusy — 1.2 x 10 4 fb, of which 71.4% is bib\ production, and 28.2% is tit\ production. The 
t\ —?• b{W branching fraction is essentially 100%, and so augments the b\ production rates. 
The &i is a quasi-stable colored particle, and will hadronize to form an i?-meson or baryon 
B ± >° = b\q or bxqq', where q is usually u or d. The properties of squark hadrons have been 
reviewed in [531 E2] • The heavy B hadron will be produced and propagate through the detector 
with minimal energy loss due to hadronic interactions. It can be of charge ±1 or 0, and in fact 
as it propagates, it has charge exchange reactions with nuclei in the detector material, so that 
the path of the 5-hadron can thus be intermittently charged or neutral [Ml [65] . 

Stable squark hadrons can be detected as muon-like events, albeit with the possibility of 
intermittently appearing and disappearing tracks. Detection can be made using dE/dx mea- 
surments in the tracking system, or time-of-flight (ToF) measurements in the muon system. 
Measuring both the B hadron momentum and velocity should allow a B mass measurement 
to about 1-2 GeV [62]. The reach of LHC for quasi-stable squark hadrons has been estimated 
in Ref. [62]. Using those results, we estimate the LHC reach with y/s = 7 TeV and 1 fb -1 
of integrated luminosity to be up to m~ b _ ~ 315 GeV, which would already test point ES4! 
The LHC reach with a/s = 14 TeV and 100 fb -1 is estimated as mg ~ 800 GeV, from just 

b\b\ production alone. Additional contributions to b\ production from cascade decays would 
increase the reach, and also make clear that 6i-pair production is not the only new physics 
process occuring at the LHC. 

We note also that many cases with a ti LMP can also be generated in ESUSY. In these 
cases, with a quasi-stable top-squark hadron, the lifetime and reach discussion is qualitatively 
similar to that given above. 
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4.4 Staus as LMP: the case of heavy quasi-stable charged particles 

(ES5) 

For ESUSY benchmark point ES5, we have fi as LMP. As with point ES4, we may again invoke 
a lighter axino to escape constraints on charged stable exotics and also mixed axion/axino CDM. 
In this case, the stau decays through loop diagrams via f\ — > ra with a lifetime typically of order 
1 sec for f a ~ 10 11 GeV, and does not significantly disrupt BBN [46]. Searches for quasi-stable 
charged particles at OPAL require m fl ~ 98.5 GeV [SB] . 

The quasi-stable stau will propagate slowly though LHC detectors much like a heavy muon, 
and leave a highly-ionizing track [61]. The dE/dx, ToF and track bending measurements should 
allow for momentum and velocity, and hence mass determinations. 

For case ES5, the total SUSY production cross section is (Jsusy — 51.4 fb. Of this, just 
about a quarter comes from stau pair production. However, since all produced sparticles now 
cascade down into the stau LMP, the total f\ production is augmented, in this case, by a factor 
four. For instance, t\t\ is produced at 26.5% rate, and t\ — > Z\t at 53.3%, t\^t tZ 2 at 7.4% and 
it ^bWx at 39.3%. Then, Z x ->• t±t? at 100% branching fraction. Also, W t -> f x v T at 100%, 
and Z% — > t^t^ at 42% and v r v T + c.c at 49%, while v T — > Z\v T) followed by further Z\ decay. 
Thus, we expect the cascade decay events to each contain at least two f\ tracks, frequently 
in company of an assortment of r-jets, together with b quarks and leptons from PU-bosons. If 
the energy of the LMP (as measured from the velocity and momentum of the track) can be 
included, E™ 1SS in such SUSY events arises only from neutrinos, and so is not especially large. 

The LHC reach for quasi-stable staus has been estimated in [62] in the case where only stau 
pair production contributes. The reach for LHC at y/s = 14 TeV with 100 fb -1 of integrated 
luminosity is found to be m T1 ~ 250 GeV. This is conservative for the present case, since fi 
production will be augmented by a factor of four from production via cascade decays of heavier 
sparticles. Moreover, as in the ES4 case, detection of more complicated events other than just 
stau pair production will bring home the richness of the new physics being detected. 

5 Summary and Conclusions 

Effective SUSY has been suggested as a model for ameliorating the flavour and CP prob- 
lems that are endemic to generic SUSY models, while maintaining naturalness in the EWSB 
sector. The general idea is that SUSY states that have large couplings to the Higgs sector 
- third generation sfermions and EW gauginos - - have masses at or below the TeV scale, 
while first and second generation sparticles that directly affect the most stringently constrained 
flavour- changing processes are heavy. The typical ESUSY spectrum consists of multi-TeV first 
and second generation sfermions along with EW gauginos, Higgs bosons and third generation 
sfermions at or below the TeV scale; gluinos may be as light as 450 GeV or as heavy as several 
TeV. 

In this paper, we have examined the LHC phenomenology of ESUSY models. Toward 
this end, we set up the ESUSY model with parameters defined at the GUT scale in Sec. [TJ 
For simplicity of incorporting flavour constraints without greatly impacting LHC physics, we 
have assumed GUT scale degeneracy of the SSB mass parameters of the first two generations, 
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but allowed the corresponding parameters for the third generation and Higgs scalars to be 
independent. 

We have delineated the viable parameter space in Sec. |2j In Sec. I3.1[ we have described the 
MCMC setup that was used for finding the portion of the entire parameter space consistent 
with various low energy constraints as well as lower limits on sparticle and Higgs boson masses 
along with theory priors that were used to obtain ESUSY spectra. Our results for the favoured 
ranges of the input parameters, assuming that the LMP is a neutralino (this is true for the bulk 
of the points), are shown by the posterior probability distributions in Fig. 0] The corresponding 
distributions for sparticle masses, illustrated in Fig. [51 indeed show the expected qualitative 
features of the ESUSY spectrum: ~TeV third generation and EW gaugino masses, 10-20 TeV 
first/second generation squarks (and sleptons, that we did not show), and gluinos ranging from 
0.5-4 TeV. The posterior probability distributions for various low energy observables and the 
thermal neutralino relic density are shown in Fig. |6j While the former generally lie close to 
their SM values (because of the nature of the ESUSY spectra), the neutralino relic density Qh 2 
ranges from below 0.01 to 10 3 with a peak value around 1 — 10. 

We note that models with large values of Q% h 2 are phenomenologically every bit as viable 
as models with a thermal neutralino LSP provided the LMP neutralino is not the true LSP, but 
decays into the LSP before the advent of Big Bang nucleosynthesis. The axino of the axion- 
axino multiplet which is present in SUSY models with the PQWW solution to the strong CP 
problem is an excellent example of an LSP that is not the LMP, and indeed models where the 
measured CDM is comprised of axions and axinos have been studied in the literature. Moreover, 
if the axino is the LSP, models where the LMP is charged or coloured are also allowed provided 
again that the LMP decays before nucleosynthesis. 

Our main results on the LHC phenomenology of ESUSY models were presented in Sec. HI 
The most important sparticle production mechansims at the LHC are gluino pair production 
for values of mi/2 ~ 700 GeV, and third generation squark pair production. A very high 
6-jet multiplicity is the hallmark of multi-jet+multilepton+.Ej llss events within the ESUSY 
framework. Gluinos mostly decay via (real or virtual) third generation squarks (see Fig. [S]) 
to third generation quarks and EW gauginos. While chargino decay branching fractions tend 
to follow those of the jy-boson, the branching fraction for neutralinos to b quarks is often 
enhanced in these scenarios. Third generation squarks frequently also have 6-quarks as their 
decay products. Finally, ESUSY events are also often rich in t quarks whose decays contribute 
to the multi-lepton component of the signal. 

We have discussed the phenomenology of five ESUSY benchmark points introduced in 
Sec. 13.31 Of these, just the point ES3 has a thermal neutralino relic density in accord with the 
measured cold DM relic density, although here the Z\ would just make up about 1/3 of the 
DM. For the other points, the LMP (which need not even be a neutralino) must decay into 
the true LSP (which might be the axino) plus SM particles. If the LMP is a neutralino, this 
decay typically occurs well after the neutralino has passed through the LHC detectors, and 
the LHC phenomenology is essentially the same as for a stable neutralino. The LMP may also 
be a coloured squark or an electrically charged stau, as illustrated by the last two benchmark 
points that we discuss. Assuming again that the axino is the LSP, the LMP is quasi-stable 
and typically traverses the detector before it decays. In these cases, the presence of a (possibly 
intermittent, in the case of squark LMP) track of a slowly moving particle (whose velocity is 
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obtained through timing information), rather than EJp lss , will be the hallmark of SUSY events. 
To conclude, effective supersymmetry with third generation squarks and EW gauginos at or 
below the TeV scale, but 10-20 TeV first/second generation sfermions ameliorates the flavour 
and CP problems that plague many SUSY models. We have examined the LHC phenomenology 
of relic-density-consistent ESUSY scenarios and shown that it may differ qualitatively from that 
in the most frequently examined mSUGRA scenario. The ESUSY picture may also be tested 
in low energy measurements. In particular, if the branching fraction for B — > rv decay, and 
especially the muon magnetic moment show significant deviation from their SM values, the 
ESUSY picture would be strongly disfavoured. 
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